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Two families of exact simple solutions of Einstein field equations for inhomogeneous stiff cosmolo- 
gies are presented. The method to obtain the solutions is based on the introduction of auxiliary 
functions in order to cast the Einstein equations in such a way that can be explicitly integrated. 
Now, despite of the equations are mathematically equivalent to the equations obtained when the 
source of matter is a scalar field, is worth to mention that the source that we consider is not a 
scalar field but a perfect fluid with the stiff equation of state. The obtained solutions are expressed 
in terms of simple functions of the used coordinates and two families of particular solutions are 
considered. The geometrical and kinematical properties of the solutions are then analyzed and the 
parameters are restricted in order to have a physically acceptable behavior. The two particular 
solutions are of the Petrov type I, the first one being regular everywhere whereas the second one 
5^ ' presents a big-bang singularity. Now, for a particular value of one of the parameters, the second 

^ |j particular solution is a vacuum solution of the Bianchi I type that reduces to the Kasner solution. 
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I. INTRODUCTION 



The main motivation for the formulation of the standardphysical cosmology, in which the universe is described as 
isotropic and completely homogeneous in all its evolution [J [3, H, Q, has been the observable fact that the present 
universe seems to be isotropic and spatially homogeneous. However, as is shown by a series of recent observations 0, 
6, 7], our universe certainly is neither exactly homogeneous nor isotropic, and also there is not sufficient reason to 
believe that the behavior of its expantion was regular at early times. Thus, in order to understand the evolution 
and large-scale structure of the universe, it is necessary to consider a more general class of cosmologies obtained by 
removing the requirement of homogeneity and isotropy. Accordingly with the above considerations, in the last two 
decades there has been an increasing interest in the study of anisotropic and inhomogeneous cosmologies, as can be 
t^J- ' see, for instance, in the references [J d, [13, EH, EH- (See also the reviews [Til fl5l fla. [TtI fl8l Il9j for a summary of 
t-H , the main work). 

Now, between the several methods used to study the spatially inhomogeneous cosmologies, the search for exact 
solutions of Einstein field equations plays a specially important role. Indeed, due to the highly nonlinear character 
i of the field equations, the knowledge of some exact solutions is crucial for the understanding of specific qualitative 
aspects that can be used as a guide for the study of more general models. Now then, the properties of the exact 
solutions must be related and must be compared with the results obtained by other different methods, such as the 
methods of dynamical systems [2(j, the methods of the theory of perturbations [2l| or the study of the structure and 
formation of singularities [22|. 

However, due to the mathematical complexity of the generic inhomogeneous models, the study of exact inhomo- 
geneous cosmologies have been limited mainly to the case of spacetimes admitting an Abelian two-parameter group 
of isometries, G2, and were initiated with the study of Gowdy of spatially compact models [23l |24| . Now, if the 
metric admits an orthogonally transitive two-parameter group of isometries, then can be written in a blocks diagonal 
form[2l|. This kind of metrics admits two spacelike commuting Killing vector fields and, as can be easily verified [25j. 
the corresponding vacuum Einstein equations leads to a complete integrable system of partial differential equations. 

On the other hand, when we consider the non vacuum Einstein equations, the mathematical complexity makes it 
very difficult to find solutions with reasonably realistic fluids. So, in almost all the works on exact inhomogeneous 
cosmologies, the solutions are obtained by taking fluids whose equation of state is of the formp = jp and, in particular, 
with 7=1. This last case, the stiff fluid state equation, leads to a complete integrable system of partial differential 
equations and was considered by Zeldovich[26j] as a good candidate in order to describe the matter content of the 
Universe in its earlier stage. 
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In agreement with the above considerations, in the present paper we present two families of simple exact inhomo- 
geneous stiff cosmologies. The method to obtain the solutions is based on the introduction of auxiliary functions in 
order to cast the Einstein equations in such a way that can be explicitly integrated and the obtained solutions can 
be expressed in terms of simple functions of the used coordinates. Now, it is worth to mention that the equations are 
written in a way that are mathematically equivalent to the equations obtained when the source of matter is a scalar 
field. However, in spite of this mathematical equivalence, the source that we are considering is a perfect fluid with 
the stiff equation of state and not a scalar field. 

The paper is organized as follows. First, in Sec. UH we present the Einstein and matter evolution equations and 
the integration procedure to obtain an explicit general solution. Then, in Sec. IIIH we consider two families of simple 
particular solutions and analyze their geometrical and kinematical properties. Finally, in Sec. HVi we summarize our 
main results. 



II. THE EINSTEIN AND MATTER EVOLUTION EQUATIONS 



In order to study inhomogeneous stiff cosmologies, we take as the starting point the metric tensor as given by the 
line element 1211 



ds* = er 2U [e^{-dt z + dr A ) + W z dx l \ + e zu dy 



(1) 



where U, 7 and W are functions of r and t only. We also consider as the matter contents a perfect fluid with the stiff 
equation of state p = p, whose energy-momentum tensor can be written as 



T a p = p(2u a up + g a p). 
With the above choices, the Einsten equations can be cast as 

-Rq/3 = 2pu a up, 

whereas the matter evolution equations can be obatined, from the conservation law 

T af3 . B = 0, 



(2) 
(3) 
(4) 



by projecting it along the temporal and spatial directions. In order to obtain the above projections, we contract the 
equation (0| with the velocity vector u a and the "spatial projection tensor" h a p = u a up + g a p, respectively. So we 
obtain 



p,f3U p + 2pu p , p = 0, 
2pu^u a , p h^+p^ a h^ a = 0, 



where we use the condition u a h a ^ = 0. 

We now impose the irrotationality condition [27 



(-$ >(U $^)l/2' 



so that the equation Q can be cast as 



which can be identically satisfied if we choose [26 



(5) 
(6) 

(7) 
(8) 
(9) 



where F is an arbitrary function of the scalar potential Now, by using ([9]), the energy-momentum tensor can be 
cast as 



T a f3 — F 



(10) 
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in such a way that the Einstein and evolution equations can be written as 

R af 3 = F$, a $ t0 , (11) 



F$>% = (12) 



where F' = |f . 

For the metric ([T]), the equation (fT2"]) reduces to 

$ rr " *,« + W-\W, r <S>,r ~ W,t$,t) + ~ *?t) = 0. (13) 

In order to solve the above equation, we assume that $ = where ^ is a new scalar potential, in such a way that 

([TBI can be written as 



2F 



w r - v4) = o. (M) 



We now choose the functional dependence of < &( - 0) in such a way that the above equation can be linearized. In order 
to do this, we take the expression in the square brackets as equal to zero, 



V + ^F- = °> (15) 



in such a way that 



ki/j = / VFd$, (16) 



where k is an arbitrary positive constant. 

Now, by using (fT6|) . is easy to see that, for any arbitrary function the energy density p, the velocity vector 

u a , the energy- momentum T a p end the Einstein system of equations can be cast as 

P = -yVwK, (17) 
+<* (18) 



i ; ; s 

r a /3 = k 2 jv.aV'.iS - T^ffa/jVvV^ j , (19) 
= *V,a^,/3- (20) 

As we can see, if we take k — 0, the above system of equations reduce to the Einstein vacuum equations. Also, it 
is worth to mention that the equations are written in a way that are mathematically equivalent to the equations 
obtained when the source of matter is a scalar field. However, in spite of this mathematical equivalence, the source 
that we are considering is a perfect fluid with the stiff equation of state and not a scalar field. 

By using the non vanishing components of the Ricci tensor for the metric (fTJ) , the Einstein and evolution equations 
can be written as the following system of partial differential equations 

W, rr ~ W.tt = 0, (21) 
(Wty r ),r - W,t), t = 0, (22) 

(WU <r ) >r - (WU,t) tt = 0, (23) 
%tW, r + 7, r ^ t = 2WU t tU >r + k 2 Wip >t il>, r + w >tr , (24) 
1,tW, t + 7#,r = W(U 2 + U 2 r ) + i [k 2 W(^ 2 + i> 2 r ) + (Wm + W, rr )) . (25) 



4 



As we can see, equation (|2"Tj) is the classical one-dimensional wave equation, whose solutions are well known. On the 
other hand, from equations (f2"2"| and (f2"3"| we can see that U(t, r) and ip(t, r) both are solutions of the same partial 
differential equation. So, in order to simplify the above system of equations, we can take U(t, r) — ip(t, r). Finally, the 
integrability conditions of the overdetermined system ([Ml) - ([23]) are equivalent to the equations (j2"Tj) - (|2"3"|) , guaranting 
so the existence of solutions. Now, as can be see from the above system of equations, the stiff fluids equations are easy 
to integrate due to the fact that, as a consequence of the stiff equation of state p — p, the equations for the metric 
functions U and W decouple from the pressure [29j | . 

In order to solve the system (|21l) - (l23t . we first consider solutions of the equation (f2"Tj) of the general form 

W(r,t) = #(r + t) +n(r-f), (26) 

where ^ and £1 are arbitrary functions. We also consider another, linearly independent, solution of (f2"Tj) written as 

V(r,t) = V(r + t)-n(r-t). (27) 

Now, by using V(r,t) and W(r, t), we define a coordinate transformation 

(t,r,x,y) ~ (V,W,x,y), (28) 

in such a way that the line element takes the form 

ds 2 = e- 2U [e 2A (-dV 2 + dW 2 ) + W 2 dx 2 } + e 2U dy 2 7 (29) 

where A is given by 

A = 7 -iln(w£-W;5). (30) 
The above coordinate transformation leads the Einstein and evolution equations to the form 

A,w=qW(tf v +tf w ), (31) 

A, v = 2qWi>, v ilj, w , (32) 



(Wtv),v - (Wif),w),w = 0, (33) 

where q = 1 + k 2 /2. As we can see, equation Q33p is equivalent to the integrability condition of the overdetermined 
system ([3Tj) - (|32l) . 

We will now consider some simple solutions of the above system. In order to do this, we seek for solutions of |3"3"1) 
of the form %p{W, V) = A(W) + B(V) and obtain 

i> = ^ (W 2 + 2V 2 ) + a 2 In W + a 3 V, (34) 

where ax, ai and a 3 are real arbitrary constants. According with this, the solution of (|3ip - (|32p is given by 

A = ^ (alV 2 + a 2 3 + 2 ai a 3 V + ai a 2 ) W 2 

2 

+ l6 W4 + qa2 lnW + 9aia 2 V 2 + 2qa 2 a 3 V. (35) 

Now, by choosing some particular solutions of (f2"6"|) - (|2T|) . we can obtain many different families of solutions for the 
full system ([2lTl - ([25]) . So, in the next section, we will present to simple families of solutions obtained by means of a 
particularly simple form of the functions ^ and f2. 



III. TWO SIMPLE FAMILIES OF PARTICULAR SOLUTIONS 



In order to obtain simple particular solutions, we will consider the two simple functions 



AM) = (36) 
/2M) - T -^, (37) 
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and will take W and £1 as denned by 



for the first family of solutions, and 



*(r,t) = /i(r,t), (38) 

n(r,t) - / 2 (r,t), (39) 

*(r,i) = /i(r,t), (40) 

n(r,t) = -/ 2 (r,t), (41) 



for the second family of solutions. 



A. The first family of solutions 

By taking the first family of solutions, we obtain for the metric functions the expressions 

W(r,t) = r, (42) 
7(r,t) = | (a^t 2 + a\ + 2a\a 3 t + a\a 2 ) r 2 
a 2 

+q-^r 4 + qa\ In r + qa\a 2 t 2 + 2qa 2 a 3 t : (43) 
lb 

U(r, t) = ^ (r 2 + 2t 2 )+ a 2 In r + a 3 t, (44) 



whereas for the fluid density we obtain the expression 

2 ; 

(2 n *\ 2 fl^/ 1 u, 2 
ait + 2a 3 t) + a 3 — — aia 2 

and for the velocity components the expressions 



''" ' ■" ~ " ' ""' "" e 2{u ^\ (45) 



fc 

tt t = — = (ait + a 3 ) , (46) 
fc /&iT «2 



(f + f). (47, 



Now, in order to obtain a physically acceptable distribution of matter, we require that p > for any value of r and t. 
So, from the expression (145|) . is easy to see that p will be no negative everywhere only if we take a\= a 2 = 0. On the 
other hand, the requirement that the velocity vector be future oriented for any value of r and t imply that we must 
take a 3 < 0. 

By imposing the above conditions, the line element ([lj take the following form 

ds 2 = e~ 2a3t le qa 3 r2 (-dt 2 + dr 2 ) + r 2 dx 2 ] + e 2a3t dy 2 , (48) 
so that the density is given by the expression 

whereas the fluid velocity takes the form 

U a = e «3(2t-a 39 r 2 )/2^ 1; ^ ^ ( 5Q ) 

guaranteing so the comoving nature of the reference frame. 
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Now, in order to see if the solution has any curvature singularity, we compute the components of the Weyl tensor 
in the natural null tetrad of the metric [25l [30] . which are given by 



*(>(*, r) = -al{2a 3 qr-q~2)e a ^ 2t ~ a ^ r2 \ 



*a(*,r) 



- 2 a 3 (2t-a 3 qr 2 ) 
> 3 > 



t (t,r) = -\a 2 {2a 3 qr + q + 2)e a ^ 2t ~ a ^ 2 \ 



(51) 
(52) 
(53) 



so that all of them are regular everywhere. Furthermore, is easy to see that the Weyl tensor is of Petrov type I. 

The kinematical quantities of the metric can be also easily computed. Thus we obtain the following expressions for 
the acceleration, the expansion and the shear of the fluid 



a a = ajqr(0, 1,0,0), 



-(|o 3 |t+o|gr/2) 



Oil 



022 



Me 



3 



2|a 3 | 



r 2 e -^-(a 3 <jr 2 +2t) 



(54) 
(55) 

(56) 
(57) 



033 = 



4|a 3 | 



r 2 e 



<-(a 3 qr 2 -6t) 



(58) 



where all the components have been computed in the natural orthonormal tetrad of the metric. As we can see the 
kinematical quantities are regular everywhere. Also is easy to see that, as 03 < 0, the components ern and 022 are 
positive, while 033 is negative. 



B. The second family of solutions 



Now, by taking the second family of solutions, we obtain for the metric functions the expressions 

W(r,t) = t, 



7(r, t) = q (a\r 2 + a\ + 2a\azr + ai<i2) 



+q-^t 4 + qa\ lnt+ qa^r 2 
16 



t z 
~2 

2qa 2 a 3 r 7 



U(r,t) 



01 



(t 2 + 2r 2 ) +a 2 lnt + a s r, 



in such a way that the fluid density is given by 



2 



2,2 



aft 



4 t 2 

whereas the velocity components are given by 



0102 — a\ (air 



2a 3 r) 



= 2(!7-7) 



(59) 

(60) 
(61) 

(62) 



(air + a 3 ) 



ait 
~2~ 



02 

t 



(63) 



(64) 



Now, as in the first family of solutions, we require that p > for any value of r and t in order to obtain a physically 
acceptable distribution of matter. From expression (|62[) is easy to see that p will be no negative everywhere only if 
we take ai = a 3 = 0, so that the expression for the density reduces to 



h 2 n 2 

K "2 ^-2fa4-q 2 +l) 



(65) 
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Now, as (a 2 — l)/(a|) < 1 < q, we have an initial singularity and then the density decreases to zero as t — > oo. On 
the other hand, the velocity vector is given by 

u a =t »a(l-9»a)(l s 0,0,0), (66) 

where, in order to have a future oriented timelike vector, we have taken a 2 < 0. Also, we can see that the spatial 
velocity is zero and thus we again have a comoving reference frame. 
The line element can be written as follows 

ds 2 = r 2a2 [t 2qa i {-dt 2 + dr 2 ) + t 2 dx 2 ] + t 2a2 dy 2 , (67) 

so that, when q = 1 (or k = 0) we have a vacuum solution of the Bianchi I type that reduces to the Kasner solution 
[3H [H[ , which can be written as [2l|, [HJ 

ds 2 = t^-^i-dt 2 + dr 2 ) + t l+d dx 2 + t l - d dy 2 , (68) 

with the Kasner parameter given by d = 1 — 2<Z2- Now, is worth to mention that another kind o inhomogeneous stiff 
cosmologies were obtained by Patel and Dadich [34[ wich also reduce to the Kasner solution. However, in contrast 
with the solution here presented, the solution of Patel and Dadich is singularity free. 

Now, in order to see if the solution has a real initial singularity, we computed the Weyl tensor in the natural null 
tetrad of the metric [25L [30j] and obtain 

*o(*,r) = l -a 2 {2a 2 -l){a 2 q-l)t- 2 ^+ 2a -- 2 , (69) 

*a(t,r) = - l -{a 2 -l)a 2 r 2 ^+ 2a -- 2 , (70) 

*4(«,r) = l -a 2 {2a 2 ~l){a 2 q~l)t- 2 ^i+ 2a -- 2 . (71) 

The scalars constructed from the Ricci and Weyl tensors diverge as t — ► 0, which corresponds to a big-bang singularity. 
Also, is easy to see that in the algebraic classification of the Reimann tensor, the metric is of Petrov type I. 

On the other hand, the kinematical quantities for this model can also be easily computed and so, by taking a 2 < 0, we 
obtain for the non- vanishing components of the acceleration, expansion and shear of the fluid the following expressions 

a a = (0,0,0,0), (72) 
a (a 2 Q - 02 + 1) 



£(a%q— 02+1) 



(73) 



(2a 2 q-2a 2 -l)t^- a ^ 
= g , (74) 

(a 2 q + 2a 2 -2)t-^+ a ^ 
V22 = jj , (75) 

{a 2 2 q - 4a 2 + lV-(«2«-3a 2 +i) 
^33 = jj , (76) 

where all the components have been computed in the natural orthonormal tetrad of the metric. It is interesting to 
see that, as the pressure gradient is zero, the acceleration is equal to zero and thus the fluid is geodesic. On the other 
hand, as a 2 < 0, all the non-vanishing components of the shear are negative. 



IV. DISCUSSION 



We presented two simple families of exact inhomogeneous stiff cosmologies. The solutions were obtained by explicitly 
integrating the Einstein and matter evolution equations by means of the introduction of an auxiliary function that 
leads to a complete integrable system of partial differential equations. Now, it is worth to mention that the equations 
were written in a way that are mathematically equivalent to the equations obtained when the source of matter is a 
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scalar field. However, despite of this mathematical equivalence, the source that we consider was not a scalar field but 
a perfect fluid with the stiff equation of state. 

A general solution was obtained that can be expressed in terms of simple functions of the used coordinates. Then, 
two families of particular solutions were considered, and their geometrical and kinematical properties were analyzed 
and the values of the parameters were restricted in order to have a physically acceptable behavior. The two particular 
solutions are of the Petrov type I, the first one being regular everywhere whereas the second one presents a big-bang 
singularity. 

Now, for a particular value of one of the parameters, the second particular solution is a vacuum solution of the 
Bianchi I type that reduces to the Kasner solution. However, although there are some other exact inhomogeneous 
stiff cosmologies that in the vacuum reduce to the Kasner solution, as is the case with the solutions obtained by Patel 
and Dadich [3J], is worth to mention that this is not a characteristic behavior of this kind of cosmologies. Indeed, an 
example of a different behavior is shown by the first particular solution considered in this paper, which in the vacuum 
is not of the Bianchi I type and do not reduces to the kasner solution. 
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